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ABSTRACT 


The present concept of automatic throttle control, as employed in 
Navy carrier-based aircraft, was investigated. The aircraft chosen for 
study was the A-/7E. The powerplant was the TF41-A-2, a turbofan engine 
with a relatively slow throttle response in the approach power paneee 

The effects of additional inputs to the approach power compensator 
were evaluated. It was shown that a considerable increase in performance 
could be achieved through the incorporation of longitudinal feedback. 
In addition, the limitations imposed on performance by large engine lags 
were found to be much less severe for systems with longitudinal feedback. 
The modifications suggested require a redesign of the approach power 


compensator system currently in use by the Navy. 
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I, INTRODUCTION 


A. APCS PERFORMANCE OBJECTIVES 

The primary purpose of an approach power compensator system (APCS) 
in Navy carrier-based aircraft is to provide the thrust necessary to 
maintain the correct airspeed during carrier approach and landing. The 
APCS can be used either in an automatic carrier landing system (ACLS) 
or in a piloted, manual approach. Although airspeed systems exist 
[Ref. 1], the systems utilized in Navy aircraft attempt to maintain a 
constant angle of attack (AOA). Supplementary inputs, such as normal 
acceleration and elevator crossfeed, are usually employed in an effort 
to maintain the "on-speed" AOA. Since the on-speed value of AOA does 
not vary from approach to approach, there is no indexing requirement. 
In an airspeed system, however, the reference airspeed is a function 
of the aircraft configuration in the landing approach. As configuration 
variables, such as weight, vary from approach to approach, so does the 
reference airspeed. This is an important factor in favor of an AOA 
System. The pilot workload which airspeed indexing would add to an 
already demanding task probably would not be tolerated, and the designer 
of an approach power compensator system should bear this in mind. 

Military specifications reflect the objective of maintaining a 
constant AOA. Reference 2 states that the output of the APCS shall be 
proportional to error in AOA, change in normal acceleration, integral 
of AOA error, and elevator position. As is noted in Ref. 1, it appears 
that the APCS design concept is overly restricted by military specifi- 
cations. nogeuer. these specifications do point out the current Navy 


APCS concept. 





B. AIRCRAFT 

The aircraft chosen for study was the A-7E, a light-attack, carrier- 
based aircraft of approximately 20,000 pounds empty gross weight. The 
A-7E was chosen because it was one of two aircraft included in a previous 
study of Navy approach power compensator problems and requirements 
[Ref. 1] and because the author has had personal experience as an A-7A 
BELOt. 

The A-7E powerplant is a turbofan engine whose throttle response 
is characterized by a relatively long and nonlinear time constant in 
the approach thrust range. The APCS is an AOA system with inputs of 
AOA, integral of AOA, normal acceleration, and unit horizontal tail 
(UHT) movement. <A detailed description of the system is contained in 


et, 3. 


jae PREVIOUS STUDIES 

Systems Technology, Incorporated conducted a study of Navy approach 
power compensator problems and requirements [Ref. 1]. The following is 
a summary of pertinent conclusions and recommendations: 


1. A fundamental conceptual difficulty with current APCS's is 
the use of AOA feedback to constrain both AOA and airspeed, as 
the two are only in phase at low frequencies. Thus, it is 
possible to control only one of these or a linear combination 
of the two with a single control input. As a consequence, APCS 
design involves a compromise between these constraints. 


2. Airframe characteristics restrict the ability to achieve 
satisfactory performance. Engine response time is a fundamental 
problem area, as the APCS cannot augment the dynamics without 
excessively overdriving the engine. It was recommended that 

the effect of engine lag on APCS performance be investigated 

in detail with a view toward specifying necessary characteristics. 





3. Gust response within the current APCS concept will, in 

general, be poor. Gust proofing against longitudinal gusts 

with an AOA system is difficult. Gust proofing against 

vertical gusts has minimal effect, as the dominant response 

is determined by the short period characteristics. 

4, It was recommended that the effects of additional 

feedbacks to the APCS be analyzed in relation to gust 

response and to optimum performance for separate con-~ 

figurations of fully automatic and manually controlled 

systems. 

The Naval Air Test Center conducted an evaluation of the APCS in 
the A-7E [Ref. 3]. The emphasis was on optimizing the performance of 
the APCS for both manual and ACLS approaches with minimal material and 
design modifications. The recommended modifications were basically 
changes in feedback gains and the incorporation of a dual time constant 
in the UHT crossfeed circuit. A shorter time constant was used for 
nose up corrections than for nose down corrections, with the net result 


that the UHT input was more effective for nose up than for nose down 


attitude changes. 


peewOB ECTIVES 

The purpose af this research was to improve the performance of the 
APCS of the A-7E aircraft as used in the ACLS. Suetewawen Lys the 
effects of additional inputs to the APCS computer were evaluated with a 
goal of improving the current concept of APCS design. In addition, the 


effect of engine response time on APCS performance was determined. 
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II, METHOD OF ANALYSIS 


A. DESCRIPTION OF THE MODEL 

The fully automatic carrier landing mode was selected for study in 
this analysis of the APCS. The manual mode with the pilot in the loop 
was not considered, Thus, the total closed loop system consisted of 
the airframe, the engine, the approach power compensator, the longitu- 
dinal automatic flight control system (AFCS), and the automatic carrier 
landing system, 

The aircraft equations of motion were linearized about the steady 
state approach conditions in accordance with standard small perturbation 
theory, as described by Etkin [Ref. 4]. Lateral-directional dynamics 
were not included, The resultant equations are shown in Figure l. Air- 
craft stability derivatives and other parameter values were taken from 
Ref. 1 and are listed in Table I. Force and moment stability deriva- 
tives are normalized with respect to mass and moment of inertia, respec- 
tively. Hence, mass and moment of inertia do not appear explicitly in 
the equations of motion, as they do in Etkin's notation. 

The A-7E engine is the TF41-A-2 turbofan. Reference 5 contains an °: 
analysis of the thrust-power lever relationship for a range of operating 
conditions. The nominal approach thrust for the operating conditions 
as shown in Table I is 3000 pounds, In this range of thrust, the engine 
response time constant is a nonlinear function of thrust. The relation- 
ship between the engine time constant, Tos and thrust, and the relation- 
ship between thrust and power lever angle (8 oa) for the given operating 
conditions are taken from Ref, 1 and are shown in Figure 2. An average 
value of 275 pounds thrust per degree SrA was used in the model. Thus, 


the variable engine lag was the only nonlinearity in the engine model. 


1 





Figure 3 is a schematic of the APCS loop closures. In small per- 
turbation theory, AOA is defined as w/Ups in radians. Thus, the inputs 
to the APCS computer are proportional plus integral w, a, (vertical 
acceleration corrected for Seeelexenteter location), and a filtered UHT 
feedforward term. 

The APCS gains are a representative set listed as 'PAX" gains in 
Ref. 1 and do not correspond exactly to those in use in the current A-/E 
APCS configuration. In addition, several simplifications were made to 
the APCS ee A complex pole at ten rad/sec in the throttle actuator 
response was neglected. This frequency is well beyond the range in 
which the APCS is effective. Mir ortle linkage hysteresis of 0.5 degrees 
was ignored. The implementation of a dual UHT input time constant, as 
described in the introduction, was not incorporated into the model. A 
gain adjust bias exists in the APCS which dienees all APCS gains by the 
same amount to compensate for ambient air temperature effects. The 
standard day value of 1.0 was AEE. KL , the vertical acceleration gain, 
has two values, the smaller value ee oe load factors in excess of 
L.l g's. Only the larger value was incorporated into the model, 

Both the eercers AFCS and the ACLS models were taken from 
Ref. 6. As above, the AFCS and ACLS gains were representative but were 
not necessarily identical to those currently in use. A complex pole 
and a first order lag at 20 rad/sec in the AFCS response was neglected. 
The AFCS output is a UHT deflection, which is a function of attitude 
error (8-6); pitch rate, and normal acceleration. The automatic carrier 
landing system was represented by the A-7E SPN-42 longitudinal control 
equation. The output, a ship-to-aircraft pitch command, is a function 


of the aircraft altitude error, Za? measured perpendicular to the ideal 


2 





glide slope. In practice, the aircraft slant range and elevation are 
measured by the SPN-42 radar, corrected for ship motion and radar 
location, transformed into cartesian coordinates, and utilized to 
command an aircraft pitch angle. 

The AFCS and ACLS longitudinal control equations are shown in 


Figure 4, APCS, AFCS, and ACLS gains are listed in Table II. 


B. DIGITAL SIMULATION 

1. Continuous System Modeling Program 

The closed loop system was simulated on the IBM 360 digital computer 
utilizing a digital simulation language, CSMP. CSMP is an acronym for 
Continuous System Modeling Program. The program is augmented by basic 
FORTRAN and provides a set of functional blocks which simulate such 
analog components as integrators, relays, and function generators. A 
detailed description of the program is contained in Ref. 7. A descrip- 


tion of CSMP functions used in the. simulation is provided in Table III. 


2. CSMP Program Components 


A sample program and output are contained in Appendix A. The 
program consists of five major sections: the aircraft equations of 
motion and engine model, the APCS model, the AFCS model, the ACLS model, 
and a gust ane model, 

a. Aircraft Equations of Motion and Engine Model 
The equations of motion as shown in Figure 1 were 


put into state format: 


{x} 


[A] {x} + [B] {R} 
(1) 


{y} = [cl] {x} + [pv] {rR} + Ce) {y} 


1 





where: 


u ; 6. 
h 
W 
{x} = {vy} =( a, {rR} =(°7 ) (2) 
s) u 
a! E 
q ce 


The A, B, C, D, and E matrices are listed in Table IV. 

The engine modeling consisted of two parts. As mentioned in 
section II.A., a relationship of 275 pounds thrust per degree SOLA was 
assumed, The engine time constant, as shown in Figure 2, was approx- 
imated by a linear function generator of fifteen unequally spaced points. 
The thrust operating range modeled was +2000 pounds about the nominal 
operating value of 3000 pounds. 

b. APCS Model 

The APCS equations were obtained from Figure 3. The 
Emeottie actuator, power lever, and engine gains were combined into a 
Single parameter. As a result, the variable labeled "PLA" in the CSMP 
program is not the actual power lever angle. Gains and other constants 
in the feedback loops were also consolidated wherever possible. AOA 
and Meer ai of AOA were kept distinct. For ease of programming, the 


6/8. transfer function was written as 


: K ' 
rs - elwe es ne K. gs = Wi, me (3) 
6. T stl s + 1/T. 6. 6. W1 


The UHT input equations in the simulation are in the form of the right 


side. of Equation (3). 
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c. AFCS Model 
The AFCS equations were obtained from the equation given 
in Figure 4 with numerical values substituted for the various gains. 
d. ACLS Model 
The altitude error, Za» was defined as the negative of h, 
as the two are measured in opposite directions perpendicular to the 
ideal glideslope. The @/Z, transfer function, as shown in Figure 4, 


was rearranged for programming as follows: 


f 





Cc RS ot : Eyl oe 
a” I ial ae 
(ese +a 
ae Sta 
— K ° ee eek 
Z 
s(stl/K) (s +7.14s+25,0) 
where 
Ken 7 Set 2 5,0) 
= cx sca 
SEI 
P. 
3.94+t.K 34 9tak 
A ee eee 
1 3.9R t_ +1.0 0 3.9R t +1.0 
il x 8 
Finally 
) aoe + Z, 8 
c _ xX Petato beac 
e Ss +b,s +b stb e HL 
where 
_ 21 7.14 = 2580 
b, 7 + 7.14 by x + 25.0 bo 7 eae 
P Pp Pp 


iS 





The SPN-42 equations in the simulation are in the form of the right 
side of Equation (5). 
e, Gust Model 
To simulate a random gust field, the output of a random 
number generator with a uniform amplitude distribution was passed . 
through a second-order filter with a break frequency of 1.0 rad/sec and 
a damping ratio of 0./07. The theoretical Bei vain of the spectral 
properties of the resultant signal is described in Appendix B. Vertical 
and horizontal RMS gust velocities of 5.0 ft/sec were used in the simu- 
lation. To insure identical gust signals for each run, a constant step 
size was used in the numerical integrations for all runs, 
f. RMS Velocity Error Criteria 
An RMS velocity error for u and w was used as a measure 
of APCS effectiveness. It was defined as; 


L 


Mae ys 2 3 
i. = qf (uc) + w(t) “Jdt (6) 
O 


where T is the run length, 
oe Die ies et Aten sor Analysis 
A fourth-order Runge-Kutta numerical integration scheme 
with a fixed step size was used in the simulation, Figure 5 is a plot 
of normal acceleration after 0.5 seconds versus step size for an initial 
altitude error. Plots of other system variables, after 0.5 seconds and 
at other values of time, show the same trend. A step size of 0.05 


seconds was chosen for the simulation. 
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C. THEORETICAL ANALYSIS 

“The theoretical analysis was primarily a root locus study of 
various APCS configurations. To reduce the enone input, multiple 
output control system to a form adaptable to the analysis, a multiple 
loop analysis technique, documented in Ref. 8, was employed. An expla- 
nation of the method, including an example which covers all applications 
used in this study, is contained in Appendix C,. Table V is a listing 
of the open loop aircraft numerators, cross product numerators, and the 


open loop denominator used in the analysis. 
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fit, FROCEUURE 


A. DIGITAL SIMULATION RUNS 

.A standard set of runs for the digital simulation was utilized for 
each configuration of the APCS, The APCS configuration shown in Figure 
3 is referred to as the standard configuration. The length of the 
simulation runs was chosen to approximate the final approach phase of 
the carrier landing. Under SPN-42 control, two of the ACLS equation 
gains, R and Ke are functions of range for a range greater than 
6000 ft. For ranges less than 6000 ft. the values of RY and K.. are 
constant at the values shown in Table II. The ACLS equation employed 
was valid for range less than 6000 ft. to the point at which deck 
motion compensation is introduced. A run time of 30 seconds was chosen 
to approximate this phase of the approach. A set of six standard 
simulation runs was chosen to simulate system response to initial high 
and low airspeed and altitude errors and to vertical and longitudinal 


gust inputs. The set of conditions for the six runs were: 


Run Chau 5 ft/sec, w(0) = -5 ft/sec 
0) = 5 ft/sec, wily) — 95> feeec 
©) eee ee Ome 
(4). z(0) = -10 ft 
(5) u, = 5 ft/sec 
(6) She = 5 ft/sec 


Initial conditions and RMS gust velocities not specifically indicated 


itm a run Are Zero, 
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B. ANALYSIS OF THE STANDARD APCS CONFIGURATION 

To apply the analysis technique described in Section II.C, the APCS 
loop closures as shown in Figure 3 were reduced to an equivalent set 
that conformed to the block diagram format shown in Figure C.1l. The 
equivalent transfer functions were a consolidation of all the elements 
in each loop closure, including the thrust to power lever terms. The 
resultant diagram is shown in Figure 6. GC. represents the proportional 


plus integral AOA feedback transfer function and is given by 


K 


ee oom er), 2 35057 
W Up T stl S T,(st1/T,) 
KO som 
=350.,/°*5/ 5 1.389 7 TK S + ean T 
O L On a 
SEs ee eee (7) 
Uy 1,8(stl/T,) (st1/T_) 


Numerically 


, Oe One 


W TS (st2.857) (st1/T.) 


G, represents the normal acceleration feedback transfer function, The 
effect of the 6.7 ft. offset of the accelerometer on APCS performance 
has been shown to be minor in the low frequency range [Ref. 1]. Thus, 
the offset was neglected in the equivalent system, simplifying the 


analysis. 


Ko *350.7 
n 
a ee arr oe es 
Z 32,2°T (st1/T,) (st1/T,) 


(8) 


Sa. 2 
T(st1) (st1/T,) 
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Ge represents the UHT crossfeed transfer function and is given by 


e 


Ke 3500 S 


ees on ee 
ae Gere eal a) 


(9) 
50, 235s 


7 T_(st0.435) (st1/T,) 


The theoretical analysis was confined to the APCS loop closures only; 
AFCS and ACLS loop closures were not considered. Justification Fomenae 
procedure is given in Section IV.B. 

The closed loop transfer functions u/6 w/8., and sae fOmmeiie 
standard APCS configuration were netwedn The engine time constant, T.> 
was kept as a parameter, and root locus plots of the zeroes of the 
closed loop denominator and a were plotted to show the effects 
of engine lag. 

Since there are two feedback paths, as shown in Figure 6, the 
closed loop: denominator is defined as A''. It was derived by applying 
equation G.5. 


a 


vr = Ww Z 
A A+ GN + GN a (10) 
a 
NT was expanded as 
a (sw-U,s8) 
Ze 0 ea 6 
NAT NT SN UpsN, (11) 


Equation 11 was substituted into Equation 10 to yield 


Te ee W Wo 8 
A A+ Cg + SC, (Nap UpN, 2 (12) 
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The expression was cleared over a common denominator, and the numerator 


expressed in two terms, one a function of Toe to give 


U + bela 
it = 
A 5(s1) (st2.857) (S#1/T_) (13) 
where 
0 
reais 1)1(c12857) A 
y = 9/44,9138°49,2419°+10.91s 46.1539e°+1.2408- 


+0.1810s-0.007699 


The zeroes of A" were plotted as a function of 1/T,, by forming an 
artificial transfer function, G, given by 
ay 
e 


G = : 





(14) 


G represents the open loop transfer function in a unity feedback, single 
loop system. This form was adaptable to a documented root locus digital 
computer program [Ref. 9]. The closed loop characteristic equation for 


the artificial system was 


A‘ = ut Iya (15) 


‘La equal to Zero repres@mts angintinree engine lag, nence, a no thrust 
condition. Thus, the root loci originate at the open loop phugoid, 
short period, and control equation roots. The zeroes of the artificial 
closure represent the roots that would be obtained for the given APCS 
feedback gains if there were no engine delay. Figure 7 is a plot of 


the root loci. Only the loci above and on the real axis are shown. An 
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engine delay of 1.8 seconds occurs at the nominal operating thrust of 
3000 pounds. This value and the values at thrust levels of +2000 pounds 
about the nominal value are indicated on Figure 7. 


§ 
std std std 


Closed loop numerator zeroes were plotted for "0 ; ‘s , and ‘sf 
. e e g 


where the subscript "std'' refers to the loop closures for the standard 


configuration. 
u u u a uw u 
N = N + GN + GN + G N 
| 6 Z, 6 AT W 6 AT 6. aah 
std (16) 
U + se NN 
s(s+1) (st2,857) (st0.435) (st1/T_) 
where 
U = s” (st1) (s#2.857) (st0.435) Ne 
e 
eee ee alleges sane “lodge: anioeee. 
Be ha etO. eyNe al hom; 
Figure 8 is a plot of the zeroes of ee /T variation. 
e 
std 
0 0 
W W 8 
Ny Pe] x” ae we + sG wn - sG_U.N +G. N 
5, 6. oAl Z AT ZO Gaeee Ss ang 
d 
st (17) 


vest Wt 
s(st1) (s+0.435) (st1/T,) 


ZZ 





U = 8(st1)(s+0.435) Ny 


e 
V = -2,155s°-511.6s'-782.78°-343.6s-68.90s-8.753 | 
Figure 9 is a plot of the zeroes of — 1/T,, variation. 
e 
std 
u u u W u Wz 
Ny = Ny 5 GON AT ay GON AT 
& & & & 
std (18) 
5 1/T,V 
s(s+2.857) (st1) (st1/T,) 
2 u 
Ue= s (st2.857) (st1)N 
& 
6 5 4 5 2 
V = 0,05453s 4+0.2789s°+1.058s +1.716s +1.141s +1.797s-0.007686 


& 
Sia 


Figure 10 is a plot of the zeroes of — 1/T., variation. 


C. INVESTIGATION OF ADDITIONAL APCS INPUTS 

In studying the effects of additional APCS feedback variables, tT. 
was fixed at the a value of 1.8 seconds. The closed loop denomina- 
tor Bad numerators evaluated at this value are shown in Table VI. The 
additional feedback variables which were considered were combinations of 
proportional, integral and derivative of longitudinal velocity, u. As 
a result, the closed loop usu, and u/6 transfer function Zeroes were 
not affected. Subsequent to each root locus analysis, a set of digital 
Simulation runs, as described in Section III.A, were completed with the 


altered configuration. A summary of the RMS velocity errors is contained 


in Table VII. 
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Longitudinal acceleration feedback was first considered, The 


equivalent transfer function for u feedback, defined as Ge eS 


350.7K* 
u T(st1/T.) 


op) 
| 


With T. equal to 1.8 seconds, this becomes 


194 ..8K* 
[T§« -s056 (19) 


t 


The closed loop denominator, defined as A''', was formed by modifying 


A'' as follows: 


Aitt = cis er su 9 
A GoM (20) 
Beare 11 is a root locus plot of the zeroes of A''' as K; was varied, 


In a like manner, longitudinal velocity feedback was considered. 


194,8K | 
u _ s+0.56 Cale) 


The closed loop denominator was again defined as A''', 


rit 11 uo, 
A Be Be (22) 
Figure 12 is a plot of the zeroes of A''' as K, was varied. 


A combination of proportional plus integral longitudinal velocity 


was considered next. 


= oa ; 
Ce soma K CK fs) (23) 


Values of = from 0.1 to 0.5 were considered. Figure 13 is a plot of the 


Henoes Of (as K, was varied with KS equal to 0.5. 
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The final configuration considered was proportional plus derivative 


of longitudinal velocity. 


Se 
on ~  g+0.56 ese (24) 
Figure 14 is a plot of the zeroes of A''' as K; was varied with values 


of equal to 0.2, 0.5, and 1.0, The phugoid roots are indicated for 
all three values of Kos The closed loop numerator for the w/6 , transfer 


function was formed. 


W W Wu 

N = N + GWU (25) 
5. Sh u 6 AT 

rev 


The subscript 'rev'' refers to revised configuration. The resultant 


moot Loci for MS equal to 0.5 are shown in Figure 15. 


D. REVISED APCS CONFIGURATION 

Based on the results of the inves eeceien of various inputs, a 
final APCS configuration was chosen, As compared to the standard 
system, this configuration had additional inputs of proportional plus 
derivative of longitudinal velocity with specific gains of 0.5 and 
Zao for AG and K‘, respectively. The revised equivalent APCS block 
diagram is shown in Figure 16. 

With the gains fixed at the values given above, root locus plots 
for the zeroes of A''! and ‘bet Eiojie 1/T variation were plotted and 


e 
rev 


are shown in Figures 17 and 18, respectively. 
Straight-line approximations to Bode plots for u/6 ., w/8 > and 


sa were constructed for comparing the standard and revised APCS 


25 





Sontigurations. The transfer functions aré listed in Tabile VIII, And 
the magnitude portion of the Bode plots are shown in Figures 19 through 


el. 
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LV.” DISCUSStONeslDeRESULIS 


A. EFFECTS OF ENGINE RESPONSE 

Strictly speaking, the incremental root locus technique showing how 
the poles and zeroes of the APCS-equipped Drrerate vary with the non- 
linear engine time constant is not valid. The use of root locus analyses 
is valid only for linear systems. However, an incremental root locus 
technique .has been used for stability analyses of systems with single-~ 
valued continuous nonlinearities [Ref. 10]. The system response is 
approximately defined by the instantaneous location of the roots, here 
a function of the instantaneous thrust. 

The initial conditions and the gust inputs used in the simulation 
were considered to be representative of those that an APCS/aircraft 
system would be expected to encounter. Visual inspection of system 
response curves did not reveal the thrust response nonlinearity; i.e., 
the curves were approximately sinusoidal in shape, indicating that the 
system performance could be evaluated by considering a fixed value of 
st equal to 1.8 seconds. The minor effect of the nonlinearity was due 
to the low variance of AT in relation to the size of the thrust 
Operating range modeled. The dependence of the response on T. was 
greater for the revised system,as can be seen by comparing Figures 7 and 
9 with Figures 17 and 18. This was a direct result of the increased 
capability to alter APCS performance with the longitudinal feedback. 

The analysis did not provide insight concerning the specific 
effects of a large engine lag on APCS performance. However, the 


limitations imposed on performance by a large engine lag were shown to 
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be less severe for a system incorporating longitudinal feedback. As 
reported in Ref. 1, the larger feedback gains or lead compensation 
necessary to improve the standard AOA system performance result in 
excesSively overdriving the engine. The longitudinal feedback provided 
an increase in performance and a more efficient use of the engine. The 
thrust response induced was directly related to the damped natural 
frequency of the phugoid mode. As the frequency was increased, as in 
Figure 12 for increasing K , peak to peak thrust values increased as 
well. However, for configurations in which the damped natural frequency 
was not appreciably increased, as in Figures 11 and 14, peak to peak 


thrust values were less than those for the standard configuration. 


B. EFFECTS OF ADDITIONAL APCS INPUTS 
As mentioned in Section ITI.B, the theoretical analysis neglected 
AFCS/ACLS loop closures. It was shown that these closures do not 
appreciably effect airspeed response for the system under study. Phugoid 
damping coefficients and damped natural frequencies were taken from 
the root locus plots of Figures 11 through 14, Measurements of the 
period and the time to one half initial amplitude of the airspeed 
responses were then made for the closed loop simulations, which included 
AFCS and ACLS closures. The damped frequencies and damping coefficients 
calculated from these responses were within a few per cent of those 
predicted by the root Psaua-eiee. Thus, it was concluded that airspeed 
response as controlled by the APCS was not significantly altered by 
AFCS and ACLS loop closures, 
1. Longitudinal Acceleration 
Only a limited increase in performance was achieved by the 


implementation of u feedback, as only a minimal increase in phugoid damp- 


N 
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De Longitudinal Velocity 


The addition of u feedback increased the absolute phugoid 
damping but provided only a minor increase in the damping coefficient 
and, therefore, no decrease in per cent overshoot. The most important 
factor was the increased damped natural frequency. It was sone ldeeed 
necessary to keep the damped frequency near the level of that of the 
standard configuration, as a much higher value would prove unacceptable 
due to the rapid oscillations in thrust which would result. 


ae Proportional Plus Integral Longitudinal Velocity 


The inclusion of an integral term had an adverse effect on 
airspeed response or an initial airspeed error. This was predicted by 
the root locus plot as shown in Figure 13. The response to a longitud- 
inal gust field was observed to be improved, as documented in Table VII. 
This shows up in the root locus plots of Figures 10 and 13 as an approx- 
imate cancellation of the phugoid poles with yt zeroes. 

4, Longitudinal Velocity Plus Acceleration 

The inclusion of a combination of proportional plus derivative 
of longitudinal velocity in the APCS provided the best increase in over= 
all system performance. For values of NS greater than 0.5, the damped 
frequency increases substantially with the feedback gain, Se Pim hope 3 
value of Ka equal to 0.5 a large increase in damping was achieved while 
the damped frequency was held nearly constant. Improvement of APCS 
performance rapidly diminished as 7 was reduced below 0.5. 

A stable first-order pole was shown to approach the origin as 
longitudinal feedback gains were increased [Figures 12 and 14], and its 
presence was evident in the digital simulations. For the configuration 


described in Section IVB4 with an initial airspeed error, a residual 
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airspeed error of opposite sign and with a value of ten per cent of the 
initial airspeed error remained and slowly decayed after the phugoid 

had essentially damped out. The effect of the pole was considered to 

be of minor significance. As shown in Figure 15, there was a correspond- 
ing zero in the w/6., numerator. Consequently, the first order pole did 


not appear in the angle of attack response. 


C. FINAL CONFIGURATION 

The choice of the final, revised configuration was based on several 
factors: phugoid damping, damped natural frequency and gust response, 
A combination of longitudinal velocity and acceleration inputs was the 
only combination able to provide a substantial increase in phugoid 
damping while maintaining a constant damped natural frequency. The 
choice of the value of Ke was somewhat arbitrary, For a value of 2,0 
the effective longitudinal velocity feedback gain was 1.0, which was of 
the same magnitude as Ky in units of volts per ft/sec; at a higher 
value of Ke the airspeed input becomes dominant over AOA, Improved 
performance at a larger value of Re was indicated by the root locus 
plot and was verified by simulation [Table VII]. 

Longitudinal gust response was substantially improved for all 
configurations with longitudinal feedback, The most significant improve- 
ments, as indicated by the RMS velocity errors, were achieved for the 
configurations in which the Anieed natural frequency was increased over 
that of the standard configuration. The ratio of u to u feedback chosen 
represented a tradeoff between gust response and frequency criteria. 

Vertical gust response as measured in simulation was not appreciably 
altered. As that was expected prior to the study (Section I.C), there 


was no theoretical analysis attempted in that area, 
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For an initial altitude error for all configurations, the induced 
airspeed error was primarily an AOA, or w(t), transient which decays in 
five seconds. The induced airspeed error, u(t), was much smaller in 
Peat ison. Thus, the addition of the longitudinal feedbacks have a 
negligible effect on responses induced by an initial altitude error. 
Based on this result, it was concluded that airspeed response at push- 
over from level flight to glideslope acquisition would be unaltered. 

The frequency response plots of Figures 19 through 21 indicate the 
frequency range which was modified by the APCS revision. The amplitude 
plots for u/6 | and u/u, indicate this range was from 0.07 to 0.5 rad/sec, 
while only the amplitudes at low frequencies were attenuated in the 
w/6 | plot. The additional favorable effect of an approximate three dB 
attenuation at the phugoid break frequency of the revised configuration 


was not indicated on the asymptotic plots. 
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V. CONCLUSIONS AND RECOMMENDATIONS 


A. APCS CONFIGURATION 

Substantial improvement in performance is possible through the use 
of longitudinal feedback. However, only an acceleration term can be 
practically obtained with the current APCS, as it would entail only the 
implementation of an accelerometer signal. Once an airspeed feedback 
is considered, the need for an airspeed reference is required, necessitat- 
ing a major APCS redesign. An airspeed reference, manually set by the 
pilot, would not be acceptable. In addition to the added pilot respon- 
sibility mentioned in Section I, the combined AOA/airspeed system 
references would require airspeed reference accuracy greater than the 
muvot could be expected to provide. 

In summary, only limited improvement in performance is possible 
within current limitations. It is recommended that the current APCS 
design concept be considered for revision, Present day aircraft, such 
as the A-7E and the F-14A/B, have precise airspeed measuring equipment 
(air data computers) and digital computers on board. The possibility 
of utilizing those devices to accurately measure and reference airspeed 
for its incorporation into the APCS should be studied. 

The possible improvement of gust response gained by the revised 
APCS configuration is partially governed by the method of sensing air- 
speed. If an inertial device is utilized, the actual airspeed, u, is 
measured, but the term aay is measured if an air data computer is 
utilized. The effect of including the gust velocity in the sensed 
airspeed signal was to degrade performance below that gained by the 


revised configrration (e.g., the RMS value of velocity error for a 
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horizontal gust input was 4.1 ft/sec. With the gust velocity included 
in the feedback signal, the error was 5.0 ft/sec, still substantially 
improved over the original value of 6.3 ft/sec). 

Angle of attack response was not degraded by the inclusion of 


longitudinal feedback but, in fact, was improved. 


B. ENGINE RESPONSE 
The effect of the nonlinear engine time constant on APCS performance 
was shown,to be minor for the range of operating conditions considered. 
The limitations imposed on APCS performance by a large engine lag 
are much less severe for a system incorporating longitudinal velocity 
and acceleration feedback than for the standard angle of attack config- 


Meeation. 


C. THEORETICAL ANALYSIS 

Restricting the multiple loop control system analysis technique to 
the APCS loop closures was —e to be valid for an evaluation of air- 
speed control. The implementation of the method provided a systematic 


approach for the study of APCS performance. 


D. GENERALIZATION OF RESULTS 

The choice of a single APCS/aircraft system for the analysis may 
restrict the validity of some of the results to that particular system. 
However, as all current Navy APCS's use the same basic AOA system, it 
is postulated that the general concept presented is applicable to other 
aircraft. Correlation of results with other systems should be attempted. 

The study was restricted to the fully automated carrier landing 

mode, The validity of the results for a pilot-controlled attitude loop 


should be investigated. 
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TABLE I 


AIRCRAFT GEOMETRY AND STABILITY DERIVATIVES 


STABILITY 
DERIVATIVES 
X ~ 0.05435 
u 
X 0.064327 
W 
Z - 0.286953 
u 
7° 0 
W 
Z ~ 0.528871 
W 
M - 0.000165 
u 
Ms - 0.000289 
W 
M - 0.007964 
W 
q 
Xe O57 32050 
e 
Z, icy 12036 
§ 
e€ 
M,  - 2.188878 
cS 
XT 0.001317 
ZAT - 0,000250 
MT 0.000004 
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APPROACH 
PARAMETERS 
S BIS, et? 
W 24,000 Lb 
SAIN) Gtr a 
yy 
MAC 10-8458 
™m 746 slugs 
h Sea Level 
6 0.002378 SHES 
ft 
Uy 218 ft/sec 
or 12 deg 
Cues 2526) (AC 
1 6.7 ft (ahead 
x 
Ofc le 





Tabi 11 


CONTROL SYSTEMS PARAMETERS 


AEGS APCS 
1.9 rad/rad Ke 143.2. volts/rad 
e 
1.0 rad/rad-sec ie 2.3 sec 
3.0 deg/g ie 2.6 volts/unit a 
K 0.55 volts/unit-sec 
O 
; I 
Ty 0.35 sec 
ACLS 
1/5.4 | aay 30.0 volts/g 
1.0 T. 1.0 sec 
0.376 
1,0 
5 
Den 
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AYA Spies EE 


CSMP FUNCTIONAL BLOCKS 


E 
Y= 1 ae ie 


0 














wae=) LNIGRL (IC, xX) 
mo) = IC 


J. 
Integrator Laplace Transform: ms 










PY ori ¥ 





REALPL (IC, P, X) 
i¢ 1 


Laplace Transform: Psa 






Order Lag 










oe e 2 
nega 2P Po a5 Po 











CMPXPL(IC, ,IC,,P, »Py+X) “Se 


io) = 1C, 
¥(0) = IC, 


Complex Pole 






Laplace Transform: 

vot | 

Z Z 

s + 2P) Pos + Py 
Y = FUNCT (X) 













Y = AFGEN(FUNCT,X) 






Arbitrary Function Gener- 
ator(Linear Interpolation) 







UNIFORM DISTRIBUTION OF 
VARIABLE Y 


Y = RNDGEN(P) 
P = ANY ODD INTEGER 












NOISE(RANDOM NUMBER) GENER- 
ATOR WITH A UNIFORM 
DISTRIBUTION 


p(Y) = Probability distri~ 
bution Function: 





a7 














TABLE IV 
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AIRCRAFT STATE EQUATIONS MATRICES 
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TABLE V 


ATRCRAFT OPEN LOOP TRANSFER FUNCTION 
NUMERATORS AND DENOMINATOR 


st + 0.973939s> + 1197797s~ + 0.109846s + 0.070689 


= 0.00025s° + 0.0003986s" - 0,0001281s + 0.00003824 


0000004078" + 0.00000422s + 0.00000320s 


0.001317s> aE 0.001195s° + 0,002435s - 0.0001321 
eS = 483.0078° = 26,38003s. = 20,1217 
0.732836s> - 0.2727048- + 40.6515s + 33.4610 
pnoncecaee + 006859768" + 0.109846s + 0.070689 
0.0191945s" at abe 2 - 0.0194913 
0.00288013s + 0.0013343 
- 0.0191945s- - 0.635365s + 0.0194913 
- 0.00060607s - 0.00085794 


- 0.00039155s~ - 0.00012806s + 0.00003824 


0.0000001131s + 0.00000320 
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TABLE VI 


AIRCRAFT TRANSFER FUNCTION NUMERATORS AND 
DENOMINATOR FOR STANDARD APCS LOOP CLOSURES 
qT. = coms Econys 


tf -_ mee rn: a eee ~ 
A 7 SED CAE BSTV THHOT Sty 4 10-303,0.2781( 0.263) 


(0.032) (-1.099)[0.346,1.44](-2.89)| 


u F. 0.7328 
‘s, SS ee aS Teas 4 [0-504,0.1821 
std . 

(0.150)[-0.176,1.38] (-0.807) (-3.41) (49.5) 
ll =i 7 | 
No.| (sti) (s¥0.435) (640.56) { (0.389,0.257] (-0.253) 
std 

(-0, 584) (-1.05) (-32.2)} 

u a Ors, , : 
‘., = S(et1) (612.857) (610,56) {[0.308,0.716] (0.00948) 
std 


- [0,0836, 1.74] (-1.58) (-3.37)} 


; | 2 2 
tat : 
Notation: [¢ , wd represents (s +200 tw ) 
(T) represents (s-T) 


positive f and negative T represent stable roots 
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TAD Ley LL 
RMS VELOCITY ERROR SUMMARY 


APCS 
Configuration 















Standard 3406 7.40 6.26 
Long Accel 
Kr=1.0 3.03 TALON 563 
=2.0 3300 O75 5a ce 
=3.0 2290 6.5/7 5603 


Long Veloc 
i =1.0 
=2.0 
=3.0 





ProptiInt u 
K (u+%/u) 
7K 


=. 10 
u 


ProptAccel u 


Co re 
K =1.0 
K*=1.0 6.76 | 4.07 
=3,0 6,18 Sess 
K =0.5 
ex o= 1.0 6.87 4.61 
=2 0% 6.55 | 4.06 
=3,0 6.35 | 3.84 
K =0.2 
4k+=1,0 ech || ye 
=3,0 6.46 | 4.33 


* Underlined values denote revised configuration. 
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APPENDIX B 


GUST MODEL ANALYSIS 


To simulate a gust input for the CSMP simulation, a random number 
generator was used, The CSMP function RNDGEN has a uniform probability 
distribution function with possible values from 0.0 to 1.0. As the 
desired mean of the gust input was 0.0, a constant value of 0.5 was 
subtracted from the RNDGEN output. The output of the random number 
generator was passed through a CSMP second-order filter with a break 
frequency of 1.0 rad/sec and a damping coefficient of 0.707 and was 
multiplied by a constant to achieve an RMS value of 5.0. The following 
is a power spectral analysis of the gust model. 

Denoting the RNDGEN output signal in the frequency domain as X(jw), 
the filter Fourier Transform as H(jw), and the resultant gust signal as 


Y(jw), the relationship between X(jw) and Y(jw) is shown in Figure B.l. 


The output power spectral 


X( jw) 





FIGURE’ B.1 GUST MODEL BLOCK DIAGRAM 


density is equal to the input power spectral density multiplied by the 


absolute value squared of the filter Fourier Transform: 
& (jw) = |aGe)| ° 6 Ge) ae 
ve” jw) | a ( 


The mean square value of the output, or average output power, can be 


determined by evaluating the autocorrelation function, oy ee at 7 =70- 
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<< 
H} 


Py) = 9 &, Ga) de 


(2) 


l sl iy tn) d 
az J Goh &,, (io) dw 
y is defined as the RMS value of y(t). 


A typical time plot of x(t) is shown in Figure B.2. 


SOC) 
O59 


~0.5 


FIGURE B. 2 RNDGEN Output vs. Time 


The random number generator output is a discrete value at each At. The 
resultant curve through the points represents the fourth-order Runge- 
Kutta approximation of the function. By approximating the curve by a 


series of steps, a rough approximation of the autocorrelation function, 


Py (T) » was determined. 


The autocorrelation function is defined by the equation, 


O.(T) = fP x(t) x(t + 7) dt (3) 


a © 


For the step series function utilized to approximate the RNDGEN output, 


each segment, At in length, is uncorrelated with every other segment; 


that is, 


Cat) = 0 [ele arene (4) 
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To determine P67) subject to the restriction of Equation (4), x(t) 


was depicted as in Figure B.3 [Ref. B.1]. By applying Equation (3) 


x(t) 
a 


[> 
Ct 


ape 
2 


a 


FIGURE B.3 Plot of x(t) to Determine Pi.) 


to Figure B.3 and requiring that Pex (7) evaluated T = 0 be equal to the 


variance, pet) was determined. 


ee 


ee)» = x (- A ) G | Se ©) 


The power spectral density can be expressed 


@ (Iw) i: fre net" alee (6) 


~~ 


By substituting Equation (5) into Equation (6) and noting that Poe (7) 


is even function, 


ee 
@ ,. (Jw) 20, ( menu + 1) coswt 

(7) 

20° 

= 5 (1 - cosuwr) 
ew ac 
By successively applying L'Hospital's Rule 

® (0) = o ‘ht (8) 

XxX ox 


i 





@ (Iw) is periodic with a period equal to 2n/At. For a small integra- 
palen step size with the filter cutoff frequency of 1.0 rad/sec, P 6 5w) 
is essentially flat in the frequency range of interest and can be 
considered a constant at the zero-frequency value shown in Equation (8). 
For a second-order filter with a break frequency of 1.0 rad/sec, the 
magnitude of the transfer function, H(jw), squared has a -40 dB/decade 
slope for frequencies greater than 1.0 rad/sec. Based on the flat power 
spectral density at low frequencies and the filter cutoff frequency, 

the evaluation of Equation (2) was approximated by taking 4. (JW) as a 


constant at its zero~-frequency value. The resulting integral was 


2 


9 W 
eee 5 s a dy (9) 
-2 (jw) + 2fu (jw) + w, 


and was of a form tabulated in Ref. B.2. Substituting for Om ands Gand 
evaluating, 


2 cae 
@) = 2 / e = 
Pa ) 0.3540 ~* At = y 


At y = 5.0 ft/sec and At = 0.05 sec, 


O. = 37.6 ft/sec (10) 


In actual simulation the value of 2, was strongly dependent on the 
RNDGEN used in the simulation, varying from 1.7 to 3.1 with a mean of 
2.4, all values falling below the predicted. RMS values of the input to 
the filter (the outputs of the RNDGEN's) were essentially constant at 
the predicted value. The underestimation of the value of oF required 
to achieve the desired value of oe was a direct result of the step 


approximation tc the fourth-order integration. The reason for the large 
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variance in output RMS values was not determined but was suspected to 
be caused by the short run time, The length of each run was 30 seconds, 
or 3600 integration points. In the final model, 0, was adjusted to 


Obtain the desired valwe of 5.0 ft/sec for a 


ow 
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APPENDIX C 


MULTIPLE LOOP CONTROL SYSTEM ANALYSIS 


The multiple loop analysis technique described below is derived 
and demonstrated for a general case in Ref. C.1. The technique is a 
procedure for expressing specific closed loop transfer functions of a 
multiple loop system in terms of the elements of the open loop, feed- 
forward and feedback matrices. The following is a two control input, 
three output example. 


Consider the system 


[A] (x} [B] {6} +  {c] (R} 


Bi? 212° 213 | ( 4r bias yo] 16462) Cia? Sy2 | 4 8S) 
25149 899% 8931) Xo(S) 1 = [Boys Boo] (S0CS)} 4 1 oq? Soa | (BoC) (1) 
eae 932° Bae Ng S77 aes Pap 231? “304 


Matrix elements are general and may contain operators. [A] is the open 
loop system matrix, {x} is a vector of variables, {6} is a vector of 
control variables, and {rR} is a vector of external disturbances. The 
system will contain one feedforward and two feedback paths to drive the 


control variable 6,, as depicted in Figure C.l 





FIGURE C.1 Control System Block Diagram 
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The control equation is 


6). = = hk - pee OF 


The (s) notation has been dropped for simplicity. Substituting Equation 


(2) into (1) and rearranging 


te. |g 


pie, 12 Pes 2 1 
“i aes 5 Xo 
me 327193928 893 238. se 
ey 11 ae 1 
=P | oo eo) | een | 1 a Ry o 
39" 7231 se 3a 


Equation (3) is represented by the matrix notation 
DAUR SEEN 16S ete Gi) LR 


where A'' denotes two loop closures and B' denotes one feedforward path. 


The closed loop system determinant can be expanded into a sum of four 


determinants 
sine Sie a3 Slo PTs 
Gos eet i (ee eres) | Cenc |b i ns 
9317 832° 833] Pole eo? oye 
(4) 
Sina na Poin SyT 121 eel 
+ =e 00 Sau? SyP2i? 82021 
S52. Sy) ay a3i° MP3? “2?3i[, 
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The determinant with a subscript of 1 is the open loop system determinant. 
Number 2 is the numerator of the open loop Xy/64 transfer function 
multiplied by h,. Number 3 is the numerator of the open loop X,/8, 
Meancfer function multiplied by ho. Number 4 is identically equal to 
zero, as the second and third columns are proportional. 

Ns denotes the numerator of the open loop X/6(s) transfer function. 


Using this notation, Equation (4) becomes 


x 


so) 3 
A'' = A +hN.* + h.N (5) 


Equation (5) can be expanded for the general case to any order and 


number of closures as 


n m 
Nae ee ns h, Nei (6) 
i=] j=l J 1 


L 


ae represents the feedback transfer function from 6. to ae 


NI 


b5 is defined as the closed loop numerator of the 


Xo? Xo 3 65 ~ rT 


X,/b5 transfer function for the above loop closures. 


From Equation (3) 


b + Y,6 h.b a + hb 


Sipe sie 


12 LL? 13 * PoPay 
Nea = | boo + ¥yPoy> Aon + hyboy> agg + Bobo, wo 
KorX&g2bo 7 Oy b eee oe a 


a9 * Y{Rage 839 7 Nene ace cay 
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Expanding Equation (7) 


b Y,b a 


Memeo? "13 Te eaas Be Mait 
BST] 
| Bert: Symes |" een aa as 
Xoo Xq2 5 7 $y b392 4399 A393 Yyb31% 839% 833 
Pio? 7a Page Fag Dios 8199 Mobi) 
+ [Dopo bybyys 894] + | bone Ange Hgby, | FO+O+0 (8) 
Bags hybs)> 433 b392 4399 hobs) 


Inspection of Equation (7) reveals that three determinants have dropped out. 


Equation (8) can be rewritten as 


x x 28 xXx.xX pee 
Ne = Ne i Ae : hyena i hyNe na” (9) 
2 2 1 94 9°4 
yo ESE Rg Uy 
where 
ae Bia? Pra» 443 
13 
N = DAD Aira 
85. 21° °22° *93 
i ligase 


Again, Equation (9) applies in general. 
In the following case the output variable is Xo » one of the feedback 
variables. 


os Eas Fes ea es an 


x 

i 2 
| =} 851° Poot Pay? 2297 MoP oy 
Xo» X32 5 7 6) 2319 Daot¥ ba), aggthgbs, 
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In the above notation 


x x x PSP? 4 
N a 2h Gee. hee (10) 


In Equation (10) another crossproduct was lost. This was a special case, 


but the notation of Equation (9) applies, noting that the lost term 
sy Siar. 
would have been hyNe - A numerator of this form makes no sense and is 
21. 


defined as identically zero. 


The following is a mumerator of a transfer function of an output 


variable due to an external disturbance. 


1 ae eld? “1a eoae 


x 
1 : 
ne] Be i a ee) 
Xo» Xa» 65 ad 5) Cay: a,oth, bs). azgthybs, 
5 ° xx pedis 
= N + hyNe et boN g” (11) 
1 1°) 154 


Equation (11) is similar to Equation (9), but the feedforarn term is 
missing. The notation of equation (11) applies to any external input 
for which there is no feedforward term, either for another control 


variable or for an actual external disturbance, as illustrated. 


es <x 
It should be noted that N : and N denote the same determinant; 
"6455 655) 
i.e., rows are not interchanged. The latter notation would arise if 
s 
N d were formed. 
2 


Xo» X32 55 7 6) 
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